Entropy for relaxation dynamics in granular media 
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We investigate the role of entropic concepts for the relaxation dynamics in granular systems. In 
these systems the existence of a geometrical frustration induces a drastic modification of the allowed 
phase space, which in its turn induces a dynamic behavior characterized by hierarchical relaxation 
phenomena with several time scales associated. In particular we show how, in the framework of a 
mean-field model introduced for the compaction phenomenon, there exists a free-energy-\i\i.e func- 
tional which decreases along the trajectories of the dynamics and which allows to account for the 
asymptotic behavior: e.g. density profile, segregation phenomena. Also we are able to perform 
the continuous limit of the above mentioned model which turns out to be a diffusive limit. In this 
framework one can single out two separate physical ingredients: the free- energy-like functional that 
defines the phase-space and the asymptotic states and a diffusion coefficient D{p) accounting for the 
velocity of approach to the asymptotic stationary states. PACS numbers: 65.50.-|-m, 45.70.-n 
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Granular media enter only partially in the framework 
of equilibrium statistical mechanics and hydrodynamics. 
Their dynamics constitutes a very complex problem of 
non-equilibrium which poses novel questions and chal- 
lenges to theorists and experimentalists |Q,|^. 

Generally speaking granular materials cannot be de- 
scribed as equilibrium systems neither from the config- 
urational point of view nor from the dynamical point 
of view. It is known in fact that these systems remain 
easily trapped in some metastable configurations which 
can last for long time intervals unless they are shaken or 
perturbed. A granular system may be in a number of 
different microscopic states at fixed macroscopic densi- 
ties, and, more in general, for a given ensemble of macro- 
scopic parameters. Many unusual properties are linked 
to this non trivial packing [|I| . The configurational space 
of these systems is very complex and presents a structure 
with several local minima. This structure induces a dy- 
namic behavior characterized by hierarchical relaxation 
phenomena with several associated time scales. A general 
mechanism bringing to the existence of such a structure 
is based on the concept of frustration, that for instance in 
granular media has a geometrical origin. The existence 
of complex geometrical interactions between the grains 
induces a rough landscape in the structure of the allow- 
able phase space and in the configurational entropy. In 
their turn these effects induce the need of complex coop- 
erative rearrangements which account for the very slow 
relaxation dynamics of these systems. At high densities 
(or very low temperatures for thermal systems) the sys- 
tem remains trapped in a local minimum and exhibits 
a non-ergodic behavior as well as very slow relaxations: 
the logarithmic compaction in granular media [^|-|7| or 
the Kohlrausch- Williams- Watts (KWW) relaxations in 
glassy systems g. 

In this paper we try to elucidate the role that the con- 
cept of entropy can play in the dynamics of these systems. 
In particular we show how, in the framework of a mean- 



field model introduced for the compaction phenomenon, 
there exists a free- energy-like functional which decreases 
along the trajectories of the dynamics and which allows 
to account for the asymptotic behavior: e.g. density pro- 
file, segregation phenomena. Furthermore the continuous 
limit of the above mentioned models allows us to com- 
ment on the relationship between entropic and dynamical 
effects in explaining the relaxation phenomena in granu- 
lar media |,|,|lO|. 

We consider a simple model which describes the evo- 
lution of a system of particles which hop on a lattice 
of fc = 0,...,N stacked planes, as introduced in M. In 
particular the system represents an ensemble of particles 
which can move up or down in a system of N layers in 
such a way that their total number is conserved. We 
ignore the correlations among particles rearrangements 
and problems related to the mechanical stability of the 
system. The master equation for the density on a generic 
plane fc, except for the k = plane, is given by: 

dtPk = (1 - Pk)D{pk)[pk-i -Pup + Pk+i ■Pdow7i)] + 

-Pfe[(l - Pk-l)D{pk-l)pdown + 

+ (1 - Pk+l)D{pk+l)Pup] (1) 

where pdown and p„p (with p^p + Pdown = 1) represents 
the probability for the particles to move downwards or 
upwards, respectively, among the different planes. With 
Pup and Pdown we can define the quantity x = Pup/ Pdown 
which quantifies the importance of gravity in the system. 
We can also associate to x a sort of temperature for the 
system given by T ^ l/log{l/x). We shall return to this 
point later on. 

D{pk) represents a sort of mobility for the particles 
given by the probability that the particle could find 
enough space to move. Apart from other effects it mainly 
takes onto account the geometrical effects of frustration, 
i.e. the fact that the packing prevents the free movement 
of the particles. 



Later on we shall show how the analysis reported here 
is very general and does not depend on the exact func- 
tional form chosen for D{p). For the sake of clarity before 
discussing the problem in its generality we shall consider 
a possible functional form for D{pk) suitable for the com- 
paction problem. In a naive way one could imagine a 
functional form like D{pk) — pk{^ — Pk') obtained by 
considering only the interactions with the nearest neigh- 
bors planes whose density is designed by pk' ■ It is easy 
to realize that such an approach does not account for the 
complexity of the problem where the packing at high den- 
sities creates long range correlations in the system, and, 
using this functional form, the equations show a trivial 
relaxation. A possible general form of D(pk), which can 
be seen as the outcome of a theory based on the existence 
of regions of cooperativity |1^| (as well as of a free- volume 
theory ||5|;|7|) for granular media, includes a term like 



D{pk) = Dq exp[-a/(l - pk) 



(2) 



The parameter a quantifies how much the shape and the 
dimensions of the particle frustrate its motion. Higher is 
a higher will be the geometrical frustration felt by the 
corresponding particle. For instance in the problems of 
parking of r-mers (i.e. segments of length r) on a line, 
under the hypothesis of an exponential distribution for 
the lengths of the empty (or filled) intervals, the prob- 
ability for each r-mer to find a sufficient space to land 
is exp(— r/(l — pj), where p is the occupation density on 
the line §. 

The question we want to address is whether there exists 
a variational principle driving the relaxation phenomena 
in this system and in general in granular media. In other 
words one could ask if, in analogy with what happens for 
a liquid system, there exists some free-e nerg y- like func- 
tional minimized (Lyapunov functional) |12(| by the dy- 
namical evolution. In a very general way it is possible 
to write explicitly a Lyapunov functional for the sys- 
tem of equations (|l|), say the functional which decreases 
monotonically along the trajectories of the motion. This 
functional can be cast in the form of a free-energy-like 
function: 



F = J2M':)kpk-S{pk 



(3) 



fc=0 



where the S{pk), the entropic-like contribution, and "f{x) 
have to be determined in a self-consistent way by impos- 
ing that F decreases for any exchange of particles be- 
tween two generic planes k and k + 1. Writing explicitly 
the expression for dF for a generic particle exchange and 
imposing dF < 0, for the particular choice of D{pk) we 
made, after some manipulations one gets for the func- 
tional (P) the expressions 



I-Pk 



1, 



S{pk) = Pklog{ -), 7(a:)=?og(-) 



Pk 



(4) 



Let us stress how S{pk) has exactly the same (critical) 
behavior, as pfc ^ 1, of the entropy of a one-dimensional 



continuous system, filled with bars of unitary length at a 
certain density p. The functional i^ is a concave function 
because d'^F/dpl > Vfc [|l|. 

As a consequence there exists a unique minimum for 
this functional and this minimum should correspond to 
the stationary state of the system. In the case of a 
monodisperse system we expect that the stationary state 
(i.e. the density profile for the system) is the one ob- 
tained, in a continuous process of shaking, after a very 
long transient. 

In the general case with arbitrary N, and for N -^ 
-l-cxD, it is possible to get the exact asymptotic station- 
ary solution for the density on each plane. We denote 
with M the total "mass" of the system, i.e. the maximal 
number of planes which can be completely filled. Using 
a standard Lagrange multiplier method, where one tries 
to find the extremum of F subject to to the constraint 
X]fc Pk — -^j the solution is given by the following im- 
plicit expression: 



where 



fipk) = fipo)x' 



m = Y^ 



(5) 



(6) 



and po is the density on the zero-th plane which is a 
complex function of the total mass of the system. In 
order to visualize the solution it is possible to extract the 
approximate explicit behaviors. In particular one gets: 



Pk 



Pk 



1 - 1/[{M - k) ■ log(l/a;)] for k « M (7) 



^(M-fc).log(l/x) f^^ ^ >> ^^ 



(8) 



The stationary solution tends thus to a step function 
9{k — M) in the limit a; — > 0. This behavior is very well 
verified in the experiments |14[| . 

Let us now ask what happens considering a bidisperse 
system, i.e. by considering two kind of particles defined 
by two different values of a (a^ = 1, a^ = a = 2 in the 
specific case). In the case the master equations for the 
densities p\ and p^. on a generic plane fc, except for the 
plane fc = 0, are given by: 

dtp):" = D^'''{pk)[p):'LiPup + P^'kliPdown)] + 

-p'k"[D^'"{Pk-l)Pdo^n + D''''{pk+l)pup] (9) 

where pk = Pk + oip\ is the total density on the /c-th 
plane and D'^'^{pk) are the probabilities for a particle s 
or h respectively, landing on the plane k, to fit the lo- 
cal geometrical environment. In this sense Z?'''*(p) are 
the analog of the mobilities for the two kind of parti- 
cles and they take into account the cooperative effects 
on the dynamics generated by the frustration. A possi- 
ble functional form for D^'^'^p) (see for a similar approach 
|-m), is given by D^ip) ^ (1 -p) exp[-l/(l -p)] and 
D^p}^ (1 - p)" cxp[-a/(l - p)] with a > 1. 



Let us now look for a Lyapunov functional equivalent 
to Eq.(|^) in this polydisperse case which now will have 
the form: 



F^Y.^^{x)npk-S{plpl)\ 



(10) 



k=t) 



where the S{p'l,pf,) is again an entropic-\ike contribu- 
tion. By repeating the calculation along the same lines 
as before one gets 

S{pIpI) = -pUlogipl) - 1) - piilogipi) 1) 

-(1 - PkKlogil - Pk) - 1) + ^05(1 - pk) (11) 

and 7(x) = log{l/x). 

Even in this case the functional F is concave and there 
exists a unique solution corresponding to the asymptotic 
stationary state for the system. Using a standard La- 
grange multiplier method one gets the result: 



6 h^{\-pkT-^e 1-"^ 

Pk = Pk- 



1 + aho{l - Pk) 



Q-l 



e ^~''t= 



Pk^ Pk- aPk 



(12) 



where ho = ^o(PojPo) i^ ^ constant which can be esti- 
mated by using the total mass M of the system: 



MPo^Po) 



Po) 



Pf>' 



(13) 



Pk is given by a non trivial function of k that can be ob- 
tained exactly but in a implicit form, p^ is a monoton- 
ically decreasing function with the maximum for fc = 0. 
Fig. |l| shows an example of the asymptotic stationary so- 
lution obtained from (p2|) with hg ~ 10^ corresponding to 

1.5-10-3. In 



0.82, po = l-- 



0.9 and ^ 

Po 



Mlog(l/x) 

particular it shows p\ and p^ as a function of pk ■ In this 
way the effect of segregation [ll^ of the "big" particles 
(with a = 2) on top of the "small" ones (with a = 1) is 
clearly visible. The introduction of two different weights 
for the two species, say p\;^ and P^^wn' gives rise to a 
rich phase space (ptpiPwp; ct) with different regions cor- 
responding to different behaviors for segregation. 

Eq.(pd]) makes evident how in the dynamics of gran- 
ular media entropic effects are far from being negligible 
and the global behavior is given by a complex interplay 
between gravitational and entropic- like effects. 

Let us now continue to push forward this point of view 
and ask whether the asymptotic stationary state corre- 
sponds to some sort of equilibrium state. At the equi- 
librium, i.e. asymptotically, one would expect that for 
any exchange of particles among different planes dF — 
for the free-energy-like functional (fiOj), or equivalently 
dS/dE = 1/T = const.. By defining the total energy 
^ = J2kLo ^Pk and the total entropy Stot = J2k=o ^iPk), 
for any exchange of particles of type 6 or s between the 
planes k and fc + 1 one has: 



dSt, 



dS{pk+i) dS{pk) 



dE 



dp\ 



fc+i 



dPk 



(14) 



It is easy to realize how in the asymptotic stationary state 
one has 



dS^. 



dE 



f of 1 / \ 

— ; — — log(—) — const. 
dE X 



(15) 



or, what is the same. 



dS" 



dE 



const, everywhere in the 



system no matter which kind of particles we use for its 
definition. It is then tempting to associate log{-) with 
the inverse of a temperature for the system. The analogy 
is made stronger by recalling that l/log{^) is the usual 
quantity associated to a temperature in Monte-Carlo dy- 
namics. From this point of view the relaxation process of 
this system would correspond to an equilibration proce- 

is made uniform everywhere. In the 



dS°, 



dE 



dure in which 

specific case of segregation one is forced to think that the 
system evolves in such a way that asymptotically all the 
particles have the same mobility, measured in this case 
in terms of the entropic change for every given displace- 
ment. 

Let us now generalize the results obtained so far. 
Eq.(P represents a particular case of a general class of 
equations that can be written as: 



dtPk = g{pk)[f{pk-i)pup + f {pk+i)Pdow7i)] - 

-f{pk)[g{pk-i)pdown + g{pk+i)pup] 



(16) 



where / and g are generic functions for which we only 
require / > 0, g > 0, df/dp > 0, dg/dp < 0. AU the 
results we have shown in the particular case (0) are valid 
under these general assumptions. In particular we can 
prove that there does exist a functional that decreases 
along the trajectories of the motion and its expression is 
given by Eq.(p|) with: 



SiPk) 



I dip) , 
I'Og-Ti—dp, 

f{p) 



7(2:) = log{ 



1 



(17) 



A deeper insight in the above mentioned phenomenol- 
ogy is obtained by considering the continuum limit for the 
model described by Eq.([l6[). More precisely we consider 
a diffusive limit that consists in scaling the space variable 
as -, the time variable as -^ and the drift term p^own— Pup 
as e. Therefore x = ek, t = e^fc/2, Pdown — Pup = e/?/2 
and we consider the evolution of u{x) = pik). 

We get the continuum limit by taking the Taylor ex- 
pansion of the right member of Eq.(|l^) around x = ke. 
For example p{k+l) = u{x+e) — u{x)+edxU+^e^dxxU+ 
0{e'). 

We get, formally, 

dru{x) = Pdxifg) + {gdxxl - fdxxg) + 0(e), (18) 
which, in the limit e ^ 0, gives 



dru{x) = Pdxifg) + [gdxxf - fdxxg)- 



(19) 



This is a non linear diffusion equation that may be con- 
veniently written in the following form 



drU = do, ( D{u)dx-^ ] , 



(20) 



where D = /g, ^ denotes the functional derivative of 
F with respect to m, i^ = J^{P ux — S{u))dx, and S" = 
log ( f ) • Notice that the functional F decreases with the 
dynamics induced by Eq. (Ey). One has, in fact, 

/r)l^ f FiF / 

dx—drU = j dx—dx f D{u)dx 

that, after an integration by parts, gives 



dF 
du 



(21) 



Therefore there exists a "free energy" -like functional, F, 
for Eq. (I20I ) which has exactly the same form of the func- 
tional defined for the discrete model (see Eq.(17)). We 
can notice that while the functional form of S and the 
value of (3 determine in a unique way the asymptotic 
state they are not sufficient to determine the dynamical 
behavior of the system. In particular in order to know it 
one should know the functional form of D{p). 

What we have discussed so far suggests the possibil- 
ity of introducing, for non-thermal systems as granular 
media, equilibrium concepts as free-energy, entropy and 
temperature. More precisely it is possible (in the case 
studied here) to predict the asymptotic state by means 
of the minimization of a suitable functional which can be 
constructed by entropic arguments. It is worth stressing 
how granular systems often exhibit memory and so the 
existence of a unique Lyapunov functional is not guaran- 
teed in general. In general one could expect that several 
Lyapunov functionals are associated to different station- 
ary states reached with different dynamical paths. 

Finally let us notice how Eq.(20) allows us to comment 
on the relationship between entropic and dynamical ef- 
fects in explaining the relaxation phenomena in granular 
media. Once the Lyapunov functional is known (and 
thus also the entropic properties of the system) it is pos- 
sible to predict the asymptotic state. However one can- 
not specify whether the asymptotic state is reached in 
a finite time unless one knows the connectivity proper- 
ties of the phase space, which in our case corresponds 
to know D{p). Acknowledgements: the authors wish 
to thank A. Coniglio, S. Krishnamurthy, H. J. Herrmann, 
M. Nicodemi and S. Roux for interesting discussions. VL 
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FIG. 1. Asymptotic density profile for the two species of 
particles in the Af-planes model described in the text. It is 
clearly visible the effect of segregation of the more frustrated 
(big) particles p* on top of the less frustrated (small) ones p^. 



